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ABSTRACT 

Using Hubble Space Telescope and ground-based U through K s photometry 
from the Great Observatories Origins Deep Survey, we measure the evolution of 
the luminosity function and luminosity density in the rest-frame optical (UBR) 
to z ~ 2, bridging the poorly explored "redshift desert" between z ~ 1 and 
z ~ 2. We also use deep near-infrared observations to measure the evolution 
in the rest-frame J-band to z ~ 1. Compared to local measurements from the 
SDSS, we find a brightening of the characteristic magnitude, M*, by ~ 2.1, 
~ 0.8 and ~ 0.7 mag between z ~ 0.1 and z ~ 1.9, in U, B, and R, respectively. 
The evolution of M* in the J-band is in the opposite sense, showing a dimming 
between redshifts z ~ 0.4 and z ~ 0.9 . This is consistent with a scenario in 
which the mean star formation rate in galaxies was higher in the past, while the 
mean stellar mass was lower, in qualitative agreement with hierarchical galaxy 
formation models. We find that the shape of the luminosity function is strongly 
dependent on spectral type and that there is strong evolution with redshift in the 
relative contribution from the different spectral types to the luminosity density. 
We find good agreement with previous measurements, supporting an increase 
in the 5-band luminosity density by a factor ~ 2 between the local value and 
z ~ 1, and little evolution between z ~ 1 and z ~ 2. We provide estimates of 
the uncertainty in our luminosity density measurements due to cosmic variance. 
We find good agreement in the luminosity function derived from an .R-selected 
and a ^-selected sample at z ~ 1, suggesting that optically selected surveys of 
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similar depth (R < 24) are not missing a significant fraction of objects at this 
redshift relative to a near-infrared-selected sample. We compare the rest-frame 
5-band luminosity functions from z 0-2 with the predictions of a semi-analytic 
hierarchical model of galaxy formation, and find qualitatively good agreement. 
In particular, the model predicts at least as many optically luminous galaxies at 
z ~ 1-2 as are implied by our observations. 

Subject headings: Galaxies: distance and redshifts - Galaxies: evolution - Galax- 
ies: fundamental parameters - Galaxies: high-redshift - Galaxies: luminosity 
function 

1. Introduction 

The luminosity function (hereafter LF) is one of the most fundamental metrics of the 
galaxy population, and is essential for characterizing statistical properties of galaxies and 
their evolution. Studying the LF as a function of cosmic epoch, galaxy type, and environment 
provides insights into the physical processes that shape galaxies. The measured parameters of 
the LF — the normalization, faint-end slope, and location of the characteristic luminosity or 
"knee" — are strong constraints on galaxy formation models and indeed are often considered 
the most basic test of a models' viability. When integrated over all luminosities, the LF 
provides a measure of the global luminosity density. The dependence of this quantity on 
wavelength and cosmic time allows us to probe the cosmic star formation and stellar mass 
assembly rate. 

To construct the LF, we require a complete and unbiased sample of galaxies to a given 
flux limit, with well-known selection functions and available redshifts. At any given redshift, 
the survey should probe a large enough volume to control cosmic variance, and contain 
enough galaxies in each magnitude bin to minimize Poisson noise. Recently, these problems 
have been largely overcome for measures of the local LF (z ~ 0) by several wide-area, multi- 
waveband surveys with follow-up spectroscopy. Recent local surveys include the 2dF Galaxy 
Redshift Survey (2dF-GRS; Colless et al. 2001; Norberg et al. 2002), the Sloan Digital Sky 
Survey (SDSS; Stoughton et al. 2002; Blanton et al. 2001, 2003) and the Two Micron All 
Sky Survey (2MASS; Jarrett et al. 2000; Kochanek et al. 2001; Cole et al. 2001). These 
studies have resolved the debate about the normalization and faint-end slope of the local 
optical LF that persisted for many years (Marzke et al. 1994; da Costa et al. 1994; Lin et 
al. 1996). 

In order to study the LF at high redshift, as is our goal here, the desire for wide area 
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coverage and large samples must be balanced against the need for deep enough photometry 
to probe faint- wards of the "knee" in the LF (L*). In addition, in order to construct the 
high redshift LF in a given fixed rest wavelength, we require photometry at wavelengths long- 
wards of the desired rest frame. Obtaining measures of the rest-frame optical LF at redshifts 
above z ~ 1 has therefore proved challenging, because of the difficulty of obtaining deep 
near-infrared (hereafter NIR) photometry over large enough fields. Obtaining spectroscopic 
redshifts for LF studies at high redshift is also extremely challenging. Most previous studies 
of the LF out to z ~ 1 based on spectroscopic surveys, such as the Canada France Redshift 
Survey (Lilly et al. 1995), CNOC2 (Lin et al. 1999) and the Caltech Faint Galaxy Redshift 
Survey (Cohen 2002) represent heroic efforts but were seriously limited by small number 
statistics and uncertainties due to cosmic variance. Two major efforts exploiting multi- 
object spectrographs on large telescopes are directed at obtaining redshifts for significant 
volumes out to z ~ 1.5 and will greatly improve this situation: the DEEP2 redshift survey 
on the Keck telescope (Davis et al. 2003) and the VIMOS-VLT Deep Survey (VVDS; Le 
Fevre et al. 2004). See also the Team Keck Treasury Redshift Survey (Wirth et al. (2004). 
However, even these surveys are optically (i?-band) selected and may suffer incompleteness 
for faint and/or red objects at z > 1. 

An alternative and powerful method of obtaining redshifts for statistically large sam- 
ples to deep flux levels is provided by the photometric redshift method, analogous to low 
resolution spectroscopy. Due to improved accuracy in photometric redshift techniques in 
recent years (A z = (\z phot — z spec |/(l + z spec )) ~ 0.1; Benitez et al. 2000; Mobasher et al. 
2004), sufficiently reliable redshifts are now available for statistical studies of galaxies. Large, 
deep mult i- waveband surveys using photometric redshifts to measure LFs to z < 1.5 have 
recently been completed by, e.g., the Las Campanas Infrared Survey (Chen et al. 2003) and 
the COMBO-17 project (Wolf et al. 2003). To higher redshifts (z < 3-5), but for relatively 
smaller areas, LFs have been presented by e.g., the Subaru Deep Survey (Kashikawa et al. 
2003), the FORS Deep Field (Gabasch et al. 2004) and by Poli et al. (2003). 

In this paper, we make use of observations obtained as part of the Great Observatories 
Origins Deep Survey (GOODS; Giavalisco et al. 2004) in the Chandra Deep Field South 
(CDF-S), also known as the GOODS-S field. We combine the results from a wider area 
(~ 1100 arcmin 2 ), optically selected (Rab < 24.5) catalog with a smaller area (~ 130 
arcmin 2 ) but very deep NIR selected (K s> ab < 23.2) catalog to measure galaxy luminosity 
functions and the luminosity density from z ~ 0.1-2. Using the optically selected sample, we 
construct these quantities from z ~ 0.1-1 for the overall galaxy population and for different 
spectral types. The NIR-selected catalog represents an unprecedented combination of area 
and depth, allowing us to probe deeper down the rest-frame optical luminosity function in 
the difficult "spectroscopic desert" regime of z ~ 1-2 than has been possible previously. We 
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construct the rest-frame optical LF to z ~ 2 using the ^-selected catalog. We also use the 
fT s -selected catalog to determine the evolution of the rest-frame J-band LF to z ~ 1. 

In the next Section, we present the multi-waveband dataset used in this study and the 
technique used to estimate photometric redshifts. The methods used to measure the LFs, 
incorporating photometric redshift errors using the redshift probability distributions, and 
our approach for estimating the uncertainty due to cosmic variance are discussed in §3. In 
§4 we present our measurements of the LF from z ~ 0.1-2 for the global population and the 
LF divided by spectral type, and compare our results with the predictions of a semi-analytic 
model of galaxy formation. We use these results to derive the integrated luminosity density 
from z ~ 0.1-2. We discuss our results and conclude in §5. 

Throughout this paper we adopt a flat cosmological constant dominated cosmology 
{Vttv = 0.7, VLm = 0.3) and a Hubble constant H =70 km s^ 1 Mpc -1 . All magnitudes are in 
the AB system (Oke & Gunn 1983). 

2. Multi- waveband data and photometric redshifts 

The photometric data for the CDF-S were obtained using ground-based optical and 
NIR observations from ESO (2.2m WFI, VLT-FORS1, NTT-SOFI and VLT-ISAAC) and 
CTIO (4m telescope) covering UBVRIJHK S , as well as space-based Hubble Space Telescope 
(HST) Advanced Camera for Surveys (ACS) observations in BViz. A summary of the 
observations and data reduction is given in Giavalisco et al. (2004). For the purpose of 
this study, two different photometric catalogs were produced. Using SExtractor (Bertin & 
Arnouts 1996), we constructed a WFI i?-band selected catalog including all ground-based 
photometry, and an ISAAC K s -band selected catalog including ISAAC JHK S and ACS 
BViz photometry. The i?-selected catalog covers ~ 1100 square arcmin while the K s - 
selected catalog covers ~ 130 square arcmin. Magnitudes used for calculating LFs are based 
on SExtractor MAG- AUTO, while colors used for determining photometric redshifts and 
spectral types are based on 3 arcsec diameter aperture magnitudes. Images were convolved 
to a common point-spread function (PSF) before measuring aperture magnitudes. The full- 
width at half maximum (FWHM) for the WFI i?-selected catalog is ~ 0.9 arcsec, while the 
ISAAC ^-selected catalog has a FWHM of ~ 0.45 arcsec. 

The apparent magnitude limits are R = 24.5 and K s = 23.2, resulting in 18,381 galaxies 
at 0.1 < z < 1.0 and 2,768 galaxies at 0.1 < z < 2.0 for the R- and K s -selected catalogs, 
respectively. The magnitude limits are sufficiently bright (~ 0.5 mag brighter than the mag- 
nitude where the number counts start to deviate from a pure power law) to allow accurate 
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photometric redshifts estimates and yet faint enough to make it possible to study the evolu- 
tion of the faint-end of the LF, as well as the integrated luminosity density, for a major part 
of the galaxy population. 

Using the available 13 (i?-selected catalog) and 7 (^-selected catalog) passbands, we 
estimate photometric redshifts for all galaxies in this survey to the magnitude limits given 
above. The photometric redshift code developed for this investigation is based on template- 
fitting method (e.g., Gwyn 1995; Mobasher et al. 1996) and includes a Bayesian prior based 
on an input LF similar to the approach used in, e.g., Kodama et al. (1999) and Dahlen et 
al. (2004). At each redshift, the absolute U-band magnitude of the object is calculated and 
compared to an input LF. The main effect of the absolute magnitude prior is to discriminate 
between cases where the chi-square fitting results in two probability peaks due to confusion 
between the Lyman- and the 4000A-break. The absolute magnitude that the object would 
have at the two redshift peaks can discriminate between the choices, i.e., absolute magnitudes 
being significantly brighter than M* is regarded as increasingly improbable. This prior is 
similar to the method of using a single bright magnitude cutoff for the absolute magnitude to 
select against wrongly assigned photometric redshifts, an approach used by e.g., COMBO-17 
(Wolf et al. 2003). 

In most cases, however, the prior does not affect the results and the photometric redshift 
is given by the best chi-square fit. Note also that we use a flat faint-end slope for the input 
prior LF. This is equivalent to not imposing any luminosity prior at all at faint magnitudes. 
This is important since we do not want to bias the slope of the LF that we later measure 
using the results from the photometric redshifts. Absorption due to intergalactic HI is 
included using the parameterization in Madau (1995). We use the four different template 
SEDs, consisting of E, Sbc, Scd, and Im, from Coleman et al. (1980), extended to UV 
and IR wavelengths by Bolzonella et al. (2000). We also include two starburst templates 
from Kinney et al. (1996, SB2 and SB3). Ten additional templates are constructed by 
interpolating between subsequent SEDs. 

Comparison between photometric and spectroscopic redshifts for ~ 400 galaxies in this 
sample results in an accuracy of A 2 ~ 0.12. After excluding a small number of outliers (~ 
3% with A z > 0.2), the accuracy becomes A z ~ 0.06. 

As part of the procedure for estimating the photometric redshift of individual galaxies 
by fitting the observed SEDs to templates, we also obtain the spectral type and the redshift 
probability distribution for each galaxy. These SED types are used in the following sections 
to study the evolution of the LFs for galaxies of different spectral types, while the redshift 
probability distribution for each galaxy is used to incorporate the errors in the photometric 
redshifts, which propagate through the estimated LFs (§3.1). 
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We divide the galaxies into three broad spectral types; early-types, late-types, and 
starbursts. Types with a best fitting SED dominated by the elliptical spectrum are defined 
as early-types, while galaxies dominated by one of the spiral galaxy spectra are called late- 
types. Starbursts are best fitted by one of the two starburst spectra. The approximate 
division between types in restframe color is that objects with B — V > 0.7 are early- types 
and objects with B — V < 0.25 are starbursts, while objects with intermediate color are 
late- types. 

Note that the spectral types, which we here divide into early-types, late-types, and 
starbursts, represent the average color of the galaxies. These do not necessarily have a 
one-to-one correlation to morphological types, e.g., morphological ellipticals, spirals, and 
irregulars. At least at moderate redshifts, however, the correlation appears to be strong. 

In Figure 1, we show the absolute magnitude vs. redshift relation for the i?-band and 
i^-band selected samples. The early-type, late-type and starburst galaxies are represented 
by red, green and blue color, respectively. Absolute magnitudes are calculated according to 
the recipe described in the next Section. 

3. Method 

In this section we present the technique used to estimate the LF and our procedure for 
accounting for the errors in our photometric redshifts. We also describe our approach for 
estimating the uncertainty in our results due to cosmic variance. 

3.1. Deriving the Luminosity Function 

Using the information on photometric redshift and best-fitting SED template, we cal- 
culate the rest-frame absolute magnitude in filter Y, My, from the observed apparent mag- 
nitude in filter X, mx, using the general equation 

M Y = m x - 5 log(£> L (z)/10 pc) - K XY (z, T) (1) 

where Kxy{z,T) is the K-correction at redshift z for template type T, correcting from 
observed filter X, to rest-frame filter Y, and D L (z) is the luminosity distance. A detailed 
description of how we calculate rest-frame magnitudes and the luminosity distance is given 
in the Appendix. 
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We use the 1/V max method (Schmidt 1968) to calculate the LF according to 



(2) 



where V^(Mj) is the observable comoving volume in which a galaxy i with absolute magnitude 
Mi is detectable, considering the apparent magnitude and redshift limits of the survey. The 
sum is taken over all galaxies in the magnitude range M — A(M)/2 < M, < M • A(M)/2. 
When evaluating equation (2), we set AM = 0.5 magnitudes. The comoving volume for any 
given galaxy is given by the redshift range z i>mhi - z i>max , where z^ mhi = Max(£ low , z rn+ ) and 
-Zi.max = Min(,2 h ig h , z m J). Here z\ ow and 2 high are the lower and upper redshift limits of the 
redshift bin where the LF is determined and z m+ and z m - are the redshifts limits where a 
galaxy with absolute magnitude Mj will have an apparent magnitude within the magnitude 
limits of the survey. For a survey covering an area Afl, the volume becomes, 



where — 1 — VIm — H is the Hubble constant, c is the speed of light and 
Dl is the luminosity distance (see Appendix). The 1/V max method is does not assume any 
underlying parametric form for the LF. Instead, after determining the number of objects in 
each magnitude bin (and their errors), it is possible to fit any desired functional form to the 
data. A drawback with the 1/Vmax method is that systematic biases may be introduced by 
large scale fluctuations and clustering within the observed field (de Lapparent et al. 1989). 

Deriving the LF requires determination of the absolute magnitude of each galaxy. With 
redshifts estimated from photometric redshift methods, the relatively large errors in redshifts 
propagate to errors in the luminosity distances and K-corrections, and therefore also to the 
absolute magnitudes. 

As faint galaxies are much more abundant than bright ones, a larger fraction of intrin- 
sically faint galaxies is shifted to the bright end, than intrinsically bright galaxies are moved 
to the faint-end (Chen et al. 2003). Thus, uncertainties in the redshifts may introduce 
biases in the determination of the LF, and lead to systematic errors when deriving the LF 
parameters, e.g., in determining the characteristic luminosity and the faint-end slope. 

To incorporate the redshift uncertainties, we use for each object its corresponding red- 
shift probability distribution, Pi(z), instead of a single redshift value. The probability dis- 
tribution is derived from the chi-square of the template fitting including input priors. 
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When summing the magnitude bins in equation (2), we use for each object a redshift 
range z_ < z < z + where the lower limit, Z-, corresponds to Mj + AM and the upper limit, 
z + , corresponds to Mj — AM, as calculated using equation (1). The contribution from each 
galaxy is weighted by the probability given by 

Pi= [ + Pi(z)dz. (4) 

This is expressed as 

Alternatively, the LF can be derived using a maximum likelihood (ML) approach (Sandage 
et al. 1979). An advantage with the ML method is that it is not as sensitive to large scale 
fluctuations as the 1/V max method. On the other hand, the ML method requires fixing the 
functional form of the LF which can be misleading since this form is not in general known 
a priori. Analogous to the method described above, Chen et al. (2003) describes a method 
for incorporating photometric redshift errors when deriving the LF using the ML method. 



3.2. Cosmic variance 

An important source of uncertainty in estimates of galaxy densities and related quanti- 
ties (LF, luminosity density, etc) in deep fields is the field-to-field variation due to large scale 
structure — commonly referred to as cosmic variance. It is generally difficult to estimate 
the magnitude of this variance empirically, as this requires measures of galaxy clustering on 
scales larger than the field in question and these are generally not available. In Somerville 
et al. (2004a), we presented an approach for estimating the cosmic variance for observed 
populations based on the expectations of clustering in the Cold Dark Matter (CDM) theory 
and a simple model for galaxy bias. These results were appropriate for the variance in the 
number counts (or number densities) of galaxies discussed in that paper, but are not directly 
applicable to the variance in the luminosity density, which is of interest here. Because more 
luminous galaxies presumably occupy more massive dark matter halos, and halo clustering 
is known to be a strong function of mass in the CDM theory, we must account for this 
in estimating the cosmic variance for mass- or luminosity-weighted quantities such as the 
luminosity density. We have developed a simple model to address this problem, which we 
briefly outline here. More details and comparisons with numerical simulations will be given 
in a future work (Newman & Somerville, in preparation). 

First, we compute the fractional variance for the underlying dark matter density field 
in each redshift bin used in our analysis. We approximate the bin geometry as a rectangular 
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solid, as in Newman & Davis (2002), and use the power spectrum for a standard (n — 1) 
ACDM cosmology with the same parameters assumed throughout this paper and a nor- 
malization <7g = 0.9. The resulting values of <xdm are given in Table 1 for both the WFI 
(i?-selected) and ISAAC (^-selected) catalog geometries. 

Next, we estimate an effective rest 5-band luminosity weighted effective bias in the 
following manner. We assume that the ratio of total halo mass to total 5-band light as a 
function of halo mass is given by the functional form 



M \ (M 
— \ (M) = 0.5 f — 



w c ) + \w c 



(6) 



following van den Bosch et al. (2003), and adopt the values for the parameters (77)0' M c , 71 
and 72 from their Model A. van den Bosch et al. (2003) showed that this model reproduces 
the luminosity function and luminosity dependent clustering of 5-band selected galaxies 
observed locally by the 2dF-GRS. We then define the luminosity-weighted effective bias as 

where pi is the total luminosity density and ^-(M, z) is the dark matter halo mass function 
at redshift z. We adopt M m i n = 1O 1O M . As the light-to-mass ratio declines rapidly with 
decreasing halo mass, the results do not depend sensitively on M m i n as long as it is of 
this order or smaller. The fractional root variance in the luminosity density is then simply 

°~L = ^lum^DM- 

Note that we assume that the relationship between light and mass (77) (M) does not 
change with redshift. This is equivalent to assuming that the halo occupation function does 
not change with time, so that all redshift dependence in our model is contained in the 
changing dark matter halo mass function and clustering amplitude rather than the way that 
galaxies of a given luminosity trace the underlying dark matter halos (e.g., Moustakas & 
Somerville 2002). It has been shown that this assumption is consistent with observations 
out to redshifts z ~ 1 (Yan et al. 2003; Coil et al. 2004), and the relatively constant 
luminosity density from z ~ 1-2 that our observations imply is also consistent with this 
assumption. Luminous (L ^ L*) galaxies must be positively biased (b > 1) in the ACDM 
model considered here, so the variance for the dark matter <tdM) probably underestimates the 
luminosity- weighted variance. If there is evolution in {^-) (M), it is likely to be in the sense 
that less massive halos will contribute a higher fraction of the total luminosity at high redshift 
than locally (as is found in semi-analytic models of galaxy formation), which would tend to 
decrease our estimate of &i um . Therefore, the true variance is almost certainly bracketed by 
<7 D m and a L . We provide both values in Table 1. We find that the 1-a uncertainty in the 
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dark matter density due to cosmic variance for the WFI (~ 1100 arcmin 2 ) fields is ~ 15 
percent, while the uncertainty on the luminosity density is ~ 16-19 percent. For the smaller 
ISAAC field (~ 130 arcmin 2 ), the DM uncertainty is ~ 15-30 percent, while the luminosity 
weighted values are ~ 16-45 percent. 

We also note that the empirical relationship between luminosity and DM halo mass 
that we have adopted here is specific to the J?-band. We expect that this relationship will 
be qualitatively similar, but perhaps different in detail, in other wavebands. The redshift 
evolution may also be more pronounced in shorter wavebands (e.g. the [/-band). Therefore 
we do not present cosmic variance estimates for the other wavebands, but expect that the 
estimates given for the J?-band should be fairly representative for the other bands. 

4. Results 

4.1. Luminosity Functions 

In this investigation, we use the 1/T4iax technique to calculate the LF, which is thereafter 
fitted to a Schechter function (Schechter 1976). We calculate the rest-frame LFs in the U- and 
S-bands using both the i?-selected (WFI based) and ^-selected (ISAAC based) catalogs, 
and the rest R- and J-band using only the ^-selected catalog. We adopt two sets of redshift 
bins appropriate to the volume and depth of each catalog. The WFI-based rest-frame U-, 
B— and ISAAC-based J-band LFs are determined in three redshift intervals, 0.1 < z < 0.5, 
0.5 < z < 0.75 and 0.75 < z < 1.0. The redshift s corresponding to the volume mid-point 
of each bin are £=0.39, £=0.64 and £=0.88, respectively (hereafter we use '£ to denote this 
redshift). For the optical bands, we further derive the type-dependent LFs in each bin. For 
the J-band, we do this only for the lowest redshift bin where statistics are sufficient. Using 
the /^-selected catalog, we extend the determination of the rest-frame U-, B- and i?-band 
LF to z ~ 2, using 6 bins with approximately equal comoving volume between z — 0.1 and 
z — 2. The redshift of the volume mid-points of these bins are £=0.62, £=0.98, £=1.24, 
£=1.48, £=1.69 and £=1.90, respectively. 

4-1.1. 'Quasi-local ' LF 

In the left panels of Figures 2-4, we show the composite rest-frame U-, B— and J- 
band LFs, as well as the type-specific LFs divided into three spectral types: early-types, 
late- types, and starbursts. The best-fitting Schechter LF parameters (M*, a, and 0*) are 
listed in Tables 2 — 4. The error ellipses for the Schechter function parameters for both 
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the total LFs and the different spectral type LFs are presented alongside the LFs. These 
contours represent 68.3% and 95% confidence intervals and correspond to A% 2 = 2.3 and 
Ax 2 = 6.0 above the best-fitting x 2 , respectively. 

The shape of the LFs varies dramatically between different spectral types. While the 
late-type LF mainly follows the total LF, showing a fairly steep faint-end slope, the early- 
type LF has a Gaussian shape which is significantly flatter at faint magnitudes. There is 
also an indication of an upturn of the LF at the very faintest magnitudes. We discuss this 
below. The starburst-type population also has a fairly steep faint-end slope, with number 
densities of these galaxies at the bright-end much smaller (by an order of magnitude) than 
other types. 

The composite and type- dependent LFs in the U- and 5-bands are similar, both in terms 
of the shape and faint-end slopes. There is a difference in that the early-type population is 
more dominant in B than in U, while the opposite is true for the starburst-type population. 
This difference is more evident in the next Section where we derive the fractional contribution 
to the luminosity density from different spectral types. 

The composite rest-frame J-band LF in Figure 4, has an overall shape similar to its op- 
tical counterparts. Also, the type-specific LFs have shapes similar to the optical. Compared 
to the optical, especially the ?7-band, the starburst population is relatively fainter in J, with 
the characteristic magnitude of the starburst population being ~ 3.0 mag fainter than the 
composite LF in J. In the U- and S-band, the difference between the characteristic magni- 
tude of the starburst-type LF and the composite LF is ~ 1.8 mag and ~ 2.5, respectively. 
The early- type population is also more dominant in J— than in [/-band. This difference 
is not significant when comparing the B- and J-bands. These trends are also reflected in 
the contribution by different populations to the total luminosity density in different bands, 
presented in the next section. 

At the faint-end of the early-type LF, we note a possible up-turn of the LF in all bands 
(most evident in B). A similar result was seen in the LF obtained for early-spectral-type 
galaxies from the 2dF-GRS (Madgwick et al. 2002). This up-turn may be an indication of 
an abundant population of faint early-type galaxies that starts to dominate the early-type 
LF at faint magnitudes. This may be analogous to the dwarf elliptical population that 
starts to dominate over normal ellipticals in local cluster of galaxies and galaxy groups (e.g., 
Binggeli et al. 1988). An alternative explanation could be contamination from faint red 
M-stars, whose SEDs are similar to early-types and may cause confusion, especially at faint 
magnitudes where it is more difficult to classify stars using e.g., PSFs. We will address this 
issue in a future paper (Dahlen et al., in preparation). 
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The faint-end upturn suggests that a single Schechter function does not provide a good 
representation of the early-type LF. To investigate how the upturn affects results, we also 
fit the Schechter parameters after excluding the faint population, similar to the approach 
in Madgwick et al. (2002) and Wolf et al. (2003). Results are given in Tables 1 — 3. As 
expected, we find that the faint-end slope is less steep after excluding the upturn. The effect 
is, however, small, mainly because the relative errors at the faint-end are large, giving low 
weight to these points when deriving the parameters over the full LF. 

In case of the J-band, we note that due to the upturn, it is possible to fit the LF with 
a straight line, corresponding to a pure power law. This is why there is no upper boundary 
on the early-type characteristic magnitude in the right panel of Figure 4. 

4-1.2. Luminosity Function Evolution 

In Figures 5 — 7, we show the composite rest-frame WFI-based U-, B— and ISAAC- 
based J-band LFs in three redshift bins 0.1 < z < 0.5, 0.5 < z < 0.75 and 0.75 < z < 1. 
Because we do not reach the same depth at all redshifts, we have to be careful when studying 
changes in M* and a with redshift; i.e., in order to determine a, it is desirable to reach at 
least ~ 3 mag fainter than M*. Also, due to the coupling between M* and a in equation 
(4), the determination of M* is affected by the depth of the survey. 

To address this, we use the following two approaches. First, we fit both M* and a in 
each redshift interval, using different limiting absolute magnitudes determined from the 
completeness in each bin. The resulting best-fitting Schechter functions are shown as solid 
red lines in Figures 5—7. Second, we fix the faint-end slope to that found in the 'quasi-local' 
sample (0.1 < z < 0.5 bin), and calculate the characteristic magnitude and normalization 
in other redshift intervals. The resulting Schechter functions are shown as dotted lines in 
Figures 5—7. For comparison, we show in the medium and high redshift bins the Schechter 
fit derived in the low redshift bin as gray lines. Results from the fits are also listed in Tables 
2-4. 

The characteristic magnitudes in both the U— and B-band LFs become brighter by 
~ 0.3 mag between z=0.4 and z=0.9. For the faint-end slope and the normalization, there is 
no clear trend with redshift. In the J-band, there is an opposite trend in the characteristic 
magnitude, i.e., it becomes somewhat fainter at higher redshifts. Again, there is no clear 
trend in the other parameters over this redshift range. The faint-end slope in all three bands 
is consistent with a value a ~ —1.3 1.4. 

Due to the coupling between M* and a in the Schechter function, their evolution with 
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redshift is not independent and a trend in one parameter may be correlated with the other. 
A better way to characterize the evolution of the LF, as well as the contribution to the total 
LF from the type dependent LFs, is to integrate the LF over all magnitudes to derive the 
average luminosity density. In the next Section we investigate this further. 

To extend the study of the rest-frame U-, B- and i?-band LFs to z ~ 2, we use the 
fT s -selected sample. This is important since galaxies start to drop-out of the i?-band se- 
lected surveys at z ^ 1, where this band probes wavelengths shorter than the 4000A-break. 
Therefore, using an optically selected sample beyond this redshift only detects galaxies that 
are very blue. As an example, with our limiting magnitudes R < 24.5 and K s < 23.2, we are 
able to detect galaxies with Mb ~ — 17.4 and Mb ~ — 16.9 at z — 0.5 in the two bands, 
respectively (assuming a Sbc galaxy template). At z = 1.5, we can detect galaxies with 
M B ~ —21.9 and M B ~ — 19.7, using the B— and ^-selected catalogs respectively, clearly 
showing that we reach significantly fainter rest-frame 5-band magnitudes at high redshift 
using the if s -selection. It is therefore useful to directly compare the LF derived using exactly 
the same methods on the i?-selected and ^-selected catalogs. 

In Figures 8-10, we show the evolution of the rest-frame U-, B-, and i?-band LF in six 
roughly equal comoving volume redshift bins from z — 0.1 to z — 2. In the lowest redshift 
bin, we fit M*, a, and 0*, while for the higher redshift bins we adopt the faint-end slope from 
the lowest redshift bin and fit only M* and 0*. The resulting best-fit Schechter functions 
are shown in the Figures as solid lines and are listed in Table 5. For reference, we show in 
the five highest redshift bins the best fitting Schechter function derived in the low redshift 
bin with dashed lines. 

Examining Figure 8 and Table 5, we find that the characteristic magnitude in the U- 
band brightens by ~ 0.9 mag from f~0.6 to f~1.9, which is a similar evolution with redshift 
as found in the i?-selected sample. The ?7-band also shows a strong decline in 0* at higher 
redshift. We see no strong evolution in M* with redshift in rest-frame B- and i?-bands, but 
note a decline in <fi* at higher redshift. 

We find good agreement between the i?-selected and -K^-selected luminosity functions in 
the redshift range where they overlap. E.g., at £=0.64 (i?-selected) and i=0.62 (^-selected), 
we find M* B = -21.46 and M* B = -21.43, as well as Mj) = -20.08 and M£ = -20.28, for 
the two different samples, respectively. 

In Figure 9 we compare our results with the predictions of a semi-analytic model of 
galaxy formation that has previously been shown to produce good agreement with the optical 
luminosity function at redshift zero and the rest-frame UV LF of Lyman break galaxies at z — 
3 (Somerville & Primack 1999; Somerville et al. 2001, hereafter SPF2001). The predictions 
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are based on the 'collisional starburst model' described in SPF2001 and computed using the 
same cosmology used throughout this paper (f2 M =0.3, Q,\=0.7, H =70 km s -1 Mpc -1 ). 

We find qualitatively good agreement between the observations and the model predic- 
tions at all redshifts to the completeness limit in each bin. At magnitudes fainter than we 
can observe, there is a discrepancy between the models and the extension of the Schechter 
function fits in the sense that the model systematically produces more faint galaxies than 
the Schechter function would imply. This discrepancy increases with redshift and will affect 
the comparison between the observed and predicted evolution of the luminosity density as 
discussed below. It is not yet clear whether the models over-predict the number of faint 
galaxies at high redshifts, or whether our assumption of a fixed, non-evolving, faint-end 
slope is incorrect. This demonstrates that probing to fainter magnitudes would provide a 
strong test of the models. 

In the discussion in §5 we compare our results on the LF with results taken from the 
literature. 



The luminosity density is calculated by integrating the LF over luminosity. To correct 
for incompleteness, we approximate the LF by the estimated Schechter function parameters 
and integrate this over all luminosities. The luminosity density is then given by: 



In Figure 11, we show the evolution of the luminosity density with redshift over the 0.1 < 
z < 1 range for rest-frame U-, B— (WFI-based) and J-bands (ISAAC-based). Filled circles 
show the evolution when fitting all Schechter parameters in each bin (M*, a, and 0*) while 
open triangles show the evolution assuming a fixed faint-end slope. The latter are somewhat 
offset in the x-direction for clarity. 

We find that between i~0.4 and f~0.9, there is a mild increase in the ?7-band luminosity 
density by a factor ~ 1.3. In the B— and J-bands, there is an increase in the luminosity 
density between the low and mid redshift bins and a slight decrease at higher redshifts. The 
overall trend in the B— and J-band luminosity density is consistent with being roughly 
constant over this epoch, especially when cosmic variance is considered. We obtain similar 
results, well within the la-errors, whether we leave the faint end slope free or fix it to the 
value in the lowest redshift bin. 

In Figure 11, we also plot predictions from semi-analytical models (SPF2001) as dashed 



4.2. Luminosity Density 
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lines. The models are in broad agreement with observations predicting a stronger evolution 
in the [/-band compared to longer wavelength bands, especially the J-band. The observed 
'bump' in the luminosity densities in the mid redshift bin, which is contrary to the more 
smooth model prediction, suggests that there is an over abundance in the 0.5 < z < 0.75 bin 
compared to the 0.75 < z < 1.0 bin due to cosmic variance. An over abundance in this 
redshift range in the CDFS is also found by the VVDS (Le Fevre et al. 2004), who suggests 
that there is a wall-like pattern at z ~ 0.7. 

In Figure 12, we show the fractional contribution to the luminosity density at 0.1 < 
z < 1.0 from galaxies of different spectral types. In the U- and 5-band, we show this for all 
three redshift bins, while for the J-band we plot only the lowest redshift bin. In the higher 
redshift bins, the small numbers of galaxies in the /^-selected catalog does not allow us to 
calculate fractional contributions with sufficient accuracy. Note that the redshift evolution 
in the fractions should only be marginally affected by cosmic variance since the total number 
of galaxies do not enter explicitly here. In the lowest redshift bin, we note that a large part 
of the luminosity density comes from late-type galaxies in all bands with the main difference 
being that the starburst contribution is higher and the early-type contribution is lower in 
the ?7-band compared to the other bands. In both the B— and J-bands, about 30% of the 
light comes from early-type galaxies. 

While there is only a weak trend in the evolution of the composite luminosity density 
in the U- and 5-band (Figure 11), there is significant evolution of type-specific luminosity 
densities, as shown in Figure 12. In the [/-band, the fraction of the luminosity density 
contributed by the starburst types increases by a factor ~ 2.4 over the redshift range, while 
the fractional contribution of both late-type and early- type galaxies decreases with redshift. 
In absolute terms, the starburst luminosity density increases by a factor ~ 3 over the redshift 
interval considered. The overall increase in the [/-band luminosity density is therefore mainly 
due to the increase in the contribution by the starburst population. The same trends seen in 
the [/-band are also present in the 5-band. However, while the starburst fraction increases, 
as in the [/-band (here a factor ~ 2.5), because the overall fraction of light contributed by 
starburst types is lower, this does not affect the evolution of the total luminosity density as 
strongly as in the [/-band. 

In Figure 13, we show a compilation of results from the literature for the rest frame 
.B-band luminosity density from < z < 2, along with our results based on the .R-selected 
and /^-selected catalogs. We show error-bars corresponding to the Poisson errors only 
(black error-bars), and error-bars including cosmic variance (grey error-bars), estimated as 
described in §3.2. We find that the luminosity densities derived from the R— and K s -selected 
samples show good agreement in the redshift range where they overlap. 
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We also find consistency between our results and those from literature, supporting a 
mild increase in the S-band luminosity between the local value and redshift z ~ 1. In 
the overlapping redshift range (0.3 £j z ^,1.2), there is also a good agreement between our 
results and COMBO-17 (Wolf et al. 2003). Compared to the local value from Norberg et 
al., both this investigation and COMBO-17 find an increase in the rest-B luminosity density 
by a factor ~ 1.7 between z ~ and z ~ 0.9. Over the same redshift range, the Lilly et al. 
(1996) results, when converted to our assumed cosmology, imply a stronger increase (factor 
~ 2.5), but the difference is within the errors. Note that the original increase reported by 
Lilly et al. was steeper (up by a factor of ~ 3.7 compared to the local value from Norberg 
et al. 2002) due to their use of a cosmology with q = 0.5 and Q — 1. At high redshifts 
(z > 1.2), we find somewhat lower luminosity densities compared to Connolly et al. (1997), 
however, we are in good agreement with the data point at z ~1.7 from Dickinson et al. 
(2003). Note here that both the Connolly et al. and the Dickinson et al. results are based 
on data from the Hubble Deep Field, but Dickinson et al. include NICMOS NIR data, 
analogous to the inclusion of ISAAC NIR data in our investigation, so likely resulting in 
more reliable photometric redshift estimates. Over the redshift range 0.5 ^ z <, 2, all the 
data are consistent with a constant luminosity density in the 5-band. 

Also shown in Fig. 13 is a comparison with the semi-analytic model predictions. The 
models predict a continuous decline in the luminosity density from z ~ 2 to the present. As 
seen in Fig. 9, the models agree well with the observed luminosity functions over the magni- 
tude range directly probed by the observations. This implies that the stronger evolution in 
the luminosity density predicted by the models is entirely due to a population of galaxies too 
faint to be directly constrained by our observations. Whether or not this population indeed 
exists remains to be seen. 



5. Discussion and Conclusions 

5.1. Optical bands 

Using the GOODS dataset, including both HST ACS and ground-based photometry 
from U through K s , we have used photometric redshifts to study the evolution of the LF 
and luminosity density over the redshift range < z < 2. We have derived the rest frame 
optical (U— and S-band) LF from 0.1 < z < 1 based on an i?-selected catalog, and used a 
smaller area but very deep ^-selected catalog to derive the rest U-, B- and i?-band LF to 
z ~ 2 and the NIR (J-band) LF to z ~ 1. 

At z < 1, we detect a brightening of the characteristic magnitude with increasing redshift 
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in optical bands based on the _R-band selected sample. A brightening in the [/-band is also 
found to z ~2 using the ^-selected sample, while there is no clear trend in the B- or the 
i?-bands at z < 0.7, based on the latter sample. 

In Figure 14 we compare the characteristic magnitude derived from GOODS with re- 
sults taken from the literature. The GOODS data are shown as circles (^-selected) and 
triangles (/^-selected), with blue, green and red colors representing rest-frame U-, B- and 
i?-band, respectively. Results from literature are converted to our adopted cosmology and 
AB magnitudes when necessary. 

In the [/-band (bottom panel) there is an excellent agreement between our results and 
the results from the FORS Deep Field (Gabasch et al. 2004; open diamonds). There is a 
clear trend of a brightening of M* with redshift. Compared to the 'local' value from SDSS 
(Blanton et al. 2001), we find a brightening of the characteristic magnitude by ~2.1 mag to 
z ~ 1.9. This is consistent with recent results from VVDS (Ilbert et al. 2004), who measures 
a brightening of by 1.8-2.4 magnitudes between z = 0.05 and z = 2.0. In the S-band 
(mid panel) there is also good agreement between GOODS data and the FORS Deep Field, 
as well as with results from Poli et al. (2003; open squares). The results from COMBO-17 
(Wolf et al. 2003; open circles) are consistent with the other results, but with higher scatter 
and with the two highest redshift points significantly brighter than the other measurements. 
Note, however, that one reason for the scatter in the Wolf et al. points is due to the use of 
a non-fixed faint-end slope when determining M* in COMBO-17. The covariance between 
M* and a causes a higher scatter in M* between bins compared to the case where a is fixed 
to a common value in all bins. The brightening of Mg between the local value from SDSS 
and our measurement at z ~ 1.9 is ~ 0.8 mag. The evolution is strongest at z < 0.7. At 
higher redshifts, our results are consistent with a flat evolution. Finally, in the i?-band (top 
panel), the scatter between measurements is the largest. There is no clear trend in either 
GOODS, COMBO-17 or data from Chen et al. (2003; open triangles) of a brightening of 
M* with redshift. When comparing to the local SDSS measurement, both GOODS and 
COMBO-17 find brighter characteristic magnitudes. The evolution in the GOODS data is 
~ 0.7 mag to z ~ 1.9, with to strongest evolution at z < 0.7, similar to the other bands. 

In summary, the non-local measurements show characteristic magnitudes that are brighter 
than the low redshift ('local') values taken from SDSS (Blanton et al. 2001; filled star) and 
2dF-GRS (Madgwick et al. 2002; open star), with a more significant difference in shorter 
wave-length bands. 

The faint-end slope of the LF in all bands is consistent with a ~ —1.3 — 1.4, and we 
see no evidence for strong evolution in the faint-end slope to z ~ 1 in any band. 
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We find that the shape of the LF is strongly dependent on spectral type, consistent 
with the results from 2dF-GRS (Madgwick et al. 2002) and COMBO-17 (Wolf et al. 2003). 
In particular, we note that the starburst-type population has a characteristic luminosity 
that is significantly fainter than the composite LF, and that the early-type LF has more of 
a Gaussian shape, with a possible upturn at faint magnitudes. This upturn in the early- 
type LF has also been reported by Madgwick et al. (2002) based on the 2dF-GRS. We 
also find significant evolution with redshift in the relative contribution from different types. 
Most prominently, the contribution from starburst galaxies to the total luminosity density 
increases with redshift. This is most evident in the [/-band, where the starburst fraction 
increases by a factor ~ 3 over the redshift range measured. At the same time, we note 
a significant decrease in the contribution from early-type galaxies with increasing redshift. 
Again, these results are consistent with similar findings based on COMBO-17 reported in 
Wolf et al. (2003). 

The LFs derived from our .R-selected and deep ^-selected catalogs agree well out to 
z ~ 1, indicating that at z < 1, ^-selected surveys (R < 24) are not likely to be missing 
a large population of the objects that would be selected in the if-band at K < 22. This 
is encouraging for the large upcoming R < 24 selected spectroscopic surveys such as DEEP 
and VVDS. 

Integrating our LF results using equation (8), we obtain estimates of the luminosity 
density in the rest-frame U-, B-, and Ribands to z ~ 2 and J-band to 2 ~ 1. The U- 
band shows a rather mild increase (about a factor of 1.3) between i~0.4 to z~0.9, based on 
the -R-selected sample. At z > 1, we do not detect any clear trend using the ^-selected 
sample. The B- and Riband luminosity densities are consistent with being constant to z ~ 2. 
This evolution could be significantly underestimated, however, if the faint end slope at high 
redshift is actually steeper than we have assumed. 

We find that our luminosity function determinations agree fairly well with the predic- 
tions of a semi-analytic model of galaxy formation over the magnitude and redshift range 
probed by the observations. The shape of the predicted luminosity function deviates from a 
Schechter form, especially at high redshift, producing an excess of both L > and L < R* 
galaxies. These problems are endemic to CDM-based models and are well known. It is 
clear, though, that the models produce at least as many luminous (R > R*) galaxies at high 
redshift (1 < z < 2) as are implied by our observations, while some previous works have 
suggested that semi-analytic models may have difficulty in this regard (e.g. Somerville et 
al. 2004b; Glazebrook et al. 2004 Nature, 430, 181). Because of the combined effects of a 
decrease in the excess of bright galaxies relative to a Schechter fit, a mild decrease in R*, 
and a flattening of the faint end of the LF with time, the models predict a monotonically 
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decreasing luminosity density from z ~ 2 to the present, in contrast with the rather flat lu- 
minosity density implied by our observations. However, the models do predict a much flatter 
dependence of luminosity density on time in the NIR bands than in the optical and UV, as 
seen in the data, consistent with an overall decrease in the global "specific star formation 
rate" (star formation rate per unit stellar mass) over time. 

5.2. NIR bands 

In the J-band, we find a trend where the characteristic magnitude, M}, gets fainter 
over the range £~0.4 to £~ 0.9 by AM} ~ 0.6. If we fix the faint-end slope to the value in 
the lowest redshift bin, the evolution is less significant, AM} ~ 0.3. A mild fading of M} 
is also found by Pozzetti et al. (2003), who measures AM} ~ 0.14 between redshifts £~0.5 
to £~1.05. However, the error-bars are larger than this difference and results are therefore 
not significant. An opposite trend is found by Feulner et al. (2003), who find a brightening 
AM} 0.6 between redshifts £~0.24 to £~0.48. 

In Figure 15, we compare the J-band LF derived in this paper with the LFs from 
Pozzetti et al. (2003) and Feulner et al. (2003). We have chosen bins at similar redshifts 
for this comparison, i.e., £=0.39, £=0.50 and £=0.88 for the GOODS, Pozzetti et al. and 
Feulner et al. (2003) data, respectively. Inspecting the Figure, we find an excellent agreement 
between the surveys. Despite this, when comparing at the Schechter function parameters 
derived in the different surveys, we find that numbers differ significantly. The characteristic 
magnitude is M} =-23.68, -22.93 and -22.98, for this investigation, Feulner et al. (2003) and 
Pozzetti et al. (2003), respectively. The faint-end slope is a =-1.48, -1.00 and -1.22, and the 
normalization is <fi =0.0008, 0.0026 and 0.0020 mag^/i^Mpc -3 , for the three investigations. 
However, comparing Schechter function parameters one-to-one can be misleading since there 
is a coupling between the different parameters. As an alternative, we use equation (8) 
the derive the luminosity density in the three surveys. We find log (p 1/ )=27. 23, 27.23 and 
27.20 erg s" 1 Hz^ 1 (h 70 /Mpc 3 ) for the three surveys, respectively. This agreement is excellent 
and further stress that comparisons between only one of the Schechter parameters may be 
misleading. 

Figure 15 also illustrates that the combination of depth and wide area in GOODS 
results in significantly better statistics, especially at the faint-end, compared to the other 
investigations. 

We find a mild increase in the J-band luminosity density between £~0.4 and £~0.9, 
results are, however, also consistent with being constant. 
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At longer rest-frame wavelengths, Drory et al. (2003) and Caputi et al. (2004) report 
a mild brightening of the characteristic magnitude with redshift in the rest-frame i^-band, 
opposite to the trend found here in the J-band. Deriving the redshift evolution in the rest- 
frame .fT-band using observed .fT-band relies on proper K- correct ions. Local K-corrections in 
the NIR are well known, however, any redshift dependence on the K-corrections will affect 
results. Also, a bias in the determination of may be introduced if there is a differential 
trend in the redshift evolution in M*, i.e., the brightening gets larger at shorter wave-lengths, 
as reported by Ilbert et al. (2004) who find an evolution in M* that is strongest in the U- 
band and becomes monotonically weaker at longer rest-frame wavelengths (including B-, 
V-, R-, and /-bands). As the observed K-ban.d probes shorter rest-frame wavelengths at 
higher redshifts, the differential trend described above could mimic brightening in M£- with 
redshift. 

We also note that making the simple assumption that the evolution in M* is proportional 
to wavelength, the results from Ilbert et al. (2004) suggest that we might expect a 'turn- 
over' at approximately the NIR J-band, where M* starts to become fainter with redshift, 
consistent with what is found here in the J-band. 

Therefore, further investigations in the rest-frame NIR are needed to firmly establish 
the evolution of the characteristic luminosity. Observations in the mid-infrared with the 
Spitzer Space Telescope will here be of great importance. 

While the optical light (especially [/-band) is related to the underlying star-formation, 
the NIR (e.g., the J-band) light is more directly related to the underlying stellar mass of 
the galaxies. The opposite trends observed in the optical and NIR for the characteristic 
magnitude therefore suggest a scenario in which the star-formation rate in galaxies increases 
with redshift, while the underlying stellar mass decreases; or equivalently in which the specific 
star formation rate (star formation rate divided by stellar mass) decreases with time. These 
trends are in qualitative agreement with predictions from the hierarchical clustering scenario. 
In a separate paper, we use the SEDs constructed here to estimate the stellar mass of each 
galaxy, and to directly estimate the rate of stellar mass build-up in our sample over the 
redshift interval from z ~ 2 to the present (Somerville et al., in preparation). 

The steep faint-end slope that we obtain in the J-band in our lowest redshift bin (a = 
-1.48l{];{jjj) is consistent with the slope a = — 1.22±g;|g derived by Pozzetti et al. (2003). 
However, the faint-end slope we find seems inconsistent with the one derived for nearby 
galaxies from 2MASS. E.g., Kochanek et al. (2001) found a = -1.09 ± 0.06 in the 2MASS 
.fT-band. We would not expect such a strong dependence on redshift or such a large difference 
between the J- and fT-bands. However, Kochanek et al. include 2MASS galaxies with 
M £ M* + 3 mag when fitting the Schechter parameters, while we reach significantly deeper, 
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M < M* + 6 mag. To investigate, we recalculated our Schechter parameters using a faint 
cutoff at M = M* + 3 mag. We find a = -1.15±^, consistent with the results from 2 MASS. 
The reason for the significant increase in the steepness of the faint-end slope at the faint 
magnitudes reached in this investigation, is that we start to probe the abundant population 
faint starburst galaxies. This is evident in Figure 4, which shows that the steep starburst 
population makes the composite LF turn steep at faint magnitudes. Contrary to this, the 
shallower 2MASS do not reach this abundant population and therefore find a flatter faint-end 
slope. This illustrates how the depth of the survey, together with the covariance between 
M* and a, can significantly affect the derived Schechter function parameters. Moreover, 
another investigation based on NIR data that is deeper than 2MASS but shallower than 
ours also found a steeper faint end slope (Huang et al. 2003). 
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A. Appendix 

A.l. Calculating rest-frame absolute magnitudes 

The rest-frame absolute magnitude My in a filter Y is calculated using the general 
formula 

M Y = m x - 5 logp L (z)/10pc) - K XY (z, T), (Al) 

where is the observed apparent magnitude in filter X, D L (z) is the luminosity distance 
and K X y(z,T) is the K-correction. The luminosity distance is given by 

D ^ z ) = ^in fL sinniW 2 [ X (n M (l + z'f + Q k (l + z'f + Q^dz'} (A2) 
#o l^fcl 1 ^ Jo 

where Qk is the curvature term defined as Qk = 1 — — &a and sinn is defined as sink for 
fife > and sin for fi fe < (Misner et al. 1973). If f2 fc = then sinn and the jf^l^-terms 
are set equal to one. 

We use the formalism in Kim et al. (1996) and Hogg et al. (2002) to calculate the 
k-correction at redshift z for a galaxy template T. The generalized equation for calculating 
the k-correction is 

- 2.5/o, 10 [ (i - - /rfAoAo , A(Ao)X(Ao) /rfAeAe/r([1 + , ]Ae )F(A e ) J ' (A3) 

where A G is the wavelength in the observed frame and A e is the wavelength in the emitted 
frame. X(X) and Y(X) are the filter transmission curves (including corrections for the de- 
tector QE), f T is the SED of the template T and g\ is the SED of the "standard source" 
used to normalize the magnitudes. For AB-magnitudes the standard source is a flat (con- 
stant) spectrum in frequency space, ^ I/ (i/)=constant. In wavelength space as used here, this 
corresponds to g\(\) = c\~ 2 g v {y). 

If we know the correct redshift and the correct SED representation of the observed 
galaxy, we can calculate the exact rest frame absolute magnitude with errors only consisting 
of the observed photometric errors. However, the true spectral type of the observed galaxy 
is generally not known. Instead we represent the the galaxy's SED with the best-fitting tem- 
plate derived from the photometric redshift fitting procedure. To minimize the dependence 
on the SED, we use the observed filter that best matches the redshifted rest-frame band of 
interest, i.e., we want to use an observed filter X that minimizes 

min{\\ x - Ay x (l + z)|} (A4) 



where Ax is the effective wavelength of the observed filter and Ay is the effective wavelength 
of the rest-frame filter in which we want to calculate the absolute magnitude. In case of 
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a perfect match, the k-correction is nearly independent of the assumed SED, while the 
dependence on SED increases with the distance between observed filter and redshifted rest- 
frame filter. 

When calculating absolute magnitudes in this investigation, we use the two observed 
filters X a and X b that satisfy equation (A4) and 

\x a < Ay x (1 + z) < \ Xb . (A5) 

From the apparent magnitudes in filters X a and X b (i.e. m Xa and m Xb ), we calculate 
the corresponding absolute magnitudes M Ya and M Yb using equation (Al). Thereafter we 
interpolate to get My, 

My = M Ya x(X Xb -X Y x(l + z))/(X Xb -X Xa ) + M Yb x (Ay x(l + z)-X x J/(X Xb -X Xa ) (A6) 

In cases where we do not have two filters available that straddle the desired rest-frame wave- 
length, we use the nearest observed magnitude according to equation (A4), and thereafter 
set My to either M Ya or My, depending on which is available. 

A. 2. Errors 

When calculating the LF using the modified l/V max method (equation 5), the redshift 
probability distribution for each object corresponds to a distribution in absolute magnitudes 
representing the uncertainty in the latter. However, for many applications it is desirable to 
derive an explicit measurement of the error in the derived rest-frame absolute magnitude. 
For completeness, we here give a recipe on how to derive these errors. 

Errors in the resulting absolute magnitudes (a to t) mainly come from three sources; 1) 
photometric errors (er m ), 2) errors due to redshift uncertainty (a z ), and 3) errors due to 
uncertainty in the best-fitting template (or). Here we consider the photometric errors to 
be known and calculate the remaining two. The magnitude errors due to uncertainty in the 
redshift are largely dominated by the uncertainty in the luminosity distance. We estimate 
this error using Monte Carlo simulations where we recalculate the luminosity distance after 
adding random errors to the photometric redshift. The distribution of random errors is 
assumed to have a la dispersion corresponding to the interval containing 68% of the redshift 
probability distribution. 

The error due to uncertainty in spectral type is estimated using MC simulations where 
we vary the SED when calculating the k-correction in equation (Al). In the photometric 
redshift code, we use a number of SEDs (numbered 1, 2, 3 etc.) where the templates follow 



-24- 



an evolutionary path going from early types to late types. In the simulations we assign to 
each object with a nominal best-fitting template N, a new SED with a random type in the 
range N-l/2 to N+l/2. With this 'interpolated' type we calculate the absolute magnitude 
and estimate variations caused by the change in SEDs. The Monte Carlo simulations are 
repeated 10,000 times for each object. 

Finally, to get the total error, we add the three parts in quadrature. 



2 2,2,2 
°"tot = °m + °" 2 + a T 
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Fig. 1. — Rest-frame 5-band magnitude vs. redshift for the i?-selected sample (top panel) 
and .^-selected sample (bottom panel). Galaxies with early- type, late-type and starburst 
spectral types are shown with red, green and blue dots, respectively. 
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Fig. 2. — The left panel shows the rest-frame 'quasi-local' [/-band luminosity function in the 
redshift range 0.1< z <0.5. Black dots show the composite LF, while red, green, and blue 
dots show the type-specific LFs for early-types, late-types, and starbursts, respectively. The 
types are determined from best fitting spectral templates. The right panel shows 68.3 and 
95 % error ellipses for M*and a from Schechter function fits to different populations. 
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Fig. 3. — The rest-frame 'quasi-local' -B-band luminosity function. Symbols are the same as 
in Figure 2. 
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Fig. 4. — The rest-frame 'quasi-local' J-band luminosity function. Symbols are the same 
as in Figure 2. Note that M* for the early- type population is not constrained at bright 
magnitudes since a Schechter function does not well represent the LF (see text for details). 
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Fig. 5. — Rest-frame £7-band luminosity function for redshift intervals 0.1 < z < 0.5 (left 
panel), 0.5 < z < 0.75 (mid panel), and 0.75 < z < 1.0 (right panel). Solid red lines show 
the best fit Schechter function, while dashed (black) lines show the best-fitting Schechter 
function where the faint-end slope, a, is fixed to the value measured in the lowest redshift 
bin. For comparison, we show with gray line in the mid and high redshift bins the best-fitting 
Schechter function found in the low redshift bin. 
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Fig. 6. — Rest-frame S-band luminosity function. Symbols are the same as in Figure 5. 




Fig. 7. — Rest-frame J-band luminosity function. Symbols are the same as in Figure 5. 
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Fig. 8. — Rest-frame (7-band luminosity function, based on the i^ s -band selected catalog. 
Each redshift bin contains the same comoving volume. The solid line shows the best fit 
Schechter function where the faint-end slope, a, has been fixed to the value measured in the 
lowest redshift bin. 




Fig. 9. — Rest-frame 5-band luminosity function (for more information see Figure 8). Blue 
lines show results from a semi-analytic model (see text). 




Fig. 10. — Rest-frame R— band luminosity function (for more information see Figure 8). 
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Fig. 11. — The evolution of the luminosity density as a function of redshift for rest-frame 
U-, B-, and J-bands is shown as filled circles. The luminosity densities in the two higher 
redshift bins assuming a fixed faint-end slope are shown with triangles (circles and triangles 
overlap in the lowest redshift bin). Predictions from a semi-analytical model are shown as 
dashed lines. 
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Fig. 12. — Fractional contribution from early- type galaxies (red points), late-type galaxies 
(green points), and starbursts (blue points) to the total luminosity density in [/-band, B- 
band and J-band. For the J-band, we only show results in the low redshift bin where 
statistics are sufficient for a determination of the fractions. 
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Fig. 13. — Evolution of the rest-frame 5-band luminosity density to z ~ 2. Blue dots show 
the GOODS measurements based on the i?-selected sample, while red dots show results 
derived from the i^ s -selected sample. Black error-bars represent statistical errors only, while 
gray error-bars also include cosmic variance. Results from the literature are taken from the 
2dF-GRS (Norberg et al. 2002; open circle), Autofib Redshift Survey (Ellis et al. 1996; 
open triangles), Canada- France Redshift Survey (Lilly et al. 1996; filled circles), COMBO- 
17 (Wolf et al. 2003), the original HDF (Connolly et al. 1997; filled triangles), and the 
HDF+NICMOS (Dickinson et al. 2003; stars). Dashed line shows the prediction from a 
semi-analytical model, as described in the text. 



-37- 



-l 1 1 1 1 1 r 1 1 1 1 1 1 1 1 1 1 1 1 1- 

R — band 



o 

Lf) 

»' CM 
* CM 

2 I t 



<x> 



* * } 



I I I I I I I I I I I I I I I I I I I I 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

B— band $ ^ 

<P [] $ t ri 



O CM 
CM 

o 

Lf) 

I 

* . 



o CM 

LO CN 

J I 



6 



<> 



H 1 1 1 1 1 1 1 1 1 1 1 1 1 h 



H h 



$ 0.2 0.4 0.6 0.8 1 

U — band 



1.2 1.4 1.6 1.8 



* 4 



T 



J i I i I i I i I i I i I i I i L 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

redshift 



Fig. 14. — Comparison between characteristic magnitudes in this work and values taken 
from literature. GOODS measurements are shown as triangles (i?-band selected sample) 
and filled circles (^-selected sample), with blue, green and red colors showing rest-frame 
U-, B- and i?-band, respectively. Literature measurements are taken from SDSS (Blanton 
et al. 2001; filled stars), 2dF-GRS (Madgwick et al. 2002; open star), FORS Deep Field 
(Gabasch et al. 2004; open diamonds), COMBO-17 (Wolf et al. 2003; open circles), Poli et 
al. (2003; open squares) and Chen et al. (2003; filled squares). 
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Fig. 15. — Rest-frame J— band LFs from this paper (open circles, £=0.39), Feulner et al. 
(2003, stars, 2=0.48) and Pozzetti et al. (2003, triangles, £=0.5 ). 
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Table 1. Cosmic Variance Estimates 



redshift range 


Kom (M P C 3 ) 


O'DM 


hum 




WFI (i?-selected) catalog (area = 1100 arcmin 2 ) 


0.1 < z < 0.5 


2.07 x 10 5 


0.15 


1.08 


0.16 


0.5 < z < 0.75 


3.68 x 10 5 


0.16 


1.17 


0.19 


0.75 < z < 1.0 


5.42 x 10 5 


0.15 


1.24 


0.19 


ISAAC (^-selected) catalog 


(area = 


= 130 


arcmin 2 ) 


0.10 < z < 0.82 


8.39 x 10 4 


0.15 


1.12 


0.16 


0.82 < z < 1.12 


8.49 x 10 4 


0.21 


1.27 


0.27 


1.12 < z < 1.37 


8.65 x 10 4 


0.23 


1.37 


0.32 


1.37 < z < 1.59 


8.44 x 10 4 


0.24 


1.50 


0.35 


1.59 < z < 1.80 


8.54 x 10 4 


0.26 


1.56 


0.41 


1.80 < z < 2.00 


8.42 x 10 4 


0.27 


1.64 


0.45 


Note. — Table columns are 1) 


redshift 


range; 


2) comov- 



ing volume corresponding to redshift range; 3) variance for 
dark matter halos; 4) luminosity-weighted effective bias; 5) 
variance in luminosity density. See Section 3.2 for details. 
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Table 2. 



Best-fitting Schechter function parameters for rest frame [/-band luminosity 

functions 



Redshift range 


M* - 5 log h 70 


a 


10- 4 (Mpc//i7 )- 3 mag- 1 


log(/v) 
erg s" 1 Hz _1 (h 7 o/Mpc 3 ) 


Spectral types 


0.1 < z < 0.5 


-19 98+ 08 


1 O1+0.02 

J — 3 -0.03 


42.5l 3 3 ; 4 8 


26.38±0.03 


All 




-i9.7iis-.ij 


-0 fi7+ 06 
u -°'-0.06 


14.5+il 




Early-types 




-19.66±g;i« 


-0 64+ - 09 


15 l +2 '° 




Early-types* 




-20 11+ - 10 

zu - ii -o.io 


-, o 9 +0.03 
1 - oz -0.03 






Late- types 




1 a OO+0.16 
lo ""-0 .20 


-1 03+ ' 14 


44.7±|1 




Starbursts 


0.5 < z < 0.75 


-20.08+°;°| 


, 10 +0.03 


64.9t« 


26.51±0.03 


All 




-19.44±g;ig 


-0 17+ - 11 


2i.s±i:i 




Early-types 




— 20 oa+O-OS 


-1 07+ ' 04 


35.4 + || 




Late- types 




iy.io_ 22 


-1 9 o+0.22 
± - zo -0.20 


32.01^ 




Starbursts 


0.75 < z < 1.0 


-20 S2+ 010 


-, ,-+0.08 
1 " : ''-0.07 


38.8 + 4 8 


26.51±0.05 


All 




-19 74+ - 24 


-0 69+ - 29 
u.oy_ 26 


9.0±i;* 




Early-types 




-20.471°;^ 


-1 26+ ' 10 
1 - zo -0.09 


19.2+if 




Late- types 




-19 71 +0 ' 18 
ly - ' i -0.18 


i ,1+0.19 
1 -J i _0.17 


26.8+™ 




Starbursts 



0.1 < z < 0.5 
0.5 < z < 0.75 
0.75 < z < 1.0 



-19.! 



n+0.08 
-0.08 



-20.44+8;°! 
-20 25+ ' 04 



Fit to the same faint-end slope 

All types, a = —1.31 



, n, +0.02 

1 -' 1 -0.03 



42.5 
40.4 
43.0 



+3.4 
3.8 

+0.4 
-0.4 
-1.7 
-0.9 



26.38±0.03 
26.54±0.03 
26.49±0.05 



All 
All 
All 



Note. — (1) i?-band selected sample. 
(2) Luminosity densities are derived by integrating LFs characterized by the best-fitting Schechter function parameters. Errors 
represent statistical errors. For early-types marked with an asterisk (*), we have excluded the faint up-turn at Mjj > —16.5. 



Table 3. Best-fitting Schechter function parameters for rest frame 5-band luminosity 

functions 



Redshift range Mg — 5 log hro 


a 


<t>* log(p„) 


Spectral types 






10- 4 (Mpc/fc 7 o)- a mag- 1 erg s" 1 Hz" 1 (h 70 /Mpc 3 ) 





Rest-frame B-band 



0.1 < z < 0.5 


91 99 +o.io 

-0.06 


, o 7 +0.02 

1 -°'-o.oi 




26.73±0.04 ± 0.06 


All 




-21 16+ ' 14 


_0 74+0 05 
U - ' 4 -0.05 






Early-types 




-20 qq+ ' 14 


-0 64+ ' 07 


15.21J-J 




Early-types* 




-2l.00lS;S° 


, OC-+0.02 
1 - oo -0.02 


21.2+^f 




Late-types 




-18.721°;- 


-1 02+"' 08 
-"-■ u -0.08 


42.61^° 




Starbursts 


0.5 < z < 0.75 


-21 46+ ' 06 


-1 9 9 +0.02 
-0.03 


3L8 ±29 


26.80±0.03 ± 0.07 


All 




-21 01+° 14 


_n sq+010 


16.7li'j 




Early-types 




-2i 2q +0 ' ()8 

zl - zy -0.08 


_1 14+0.04 


23.4±2-6 




Late-types 




_1Q =1+0.16 

iy.ai_ 16 


-1 04+°- 16 


27.8tg;g 




Starbursts 


0.75 < z < 1.0 


-21 46+ ' 10 


-1 24+ 08 
L - z -0.08 


24.5+|- 7 


26.70±0.05±0.07 


All 




— 21 44+ - 26 
zi - 44 -0.28 


r\ 7 9+0.28 
u - ,z -0.26 






Early-types 




-21 15+ ' 10 


-0.87l° ; 12 


22.0l| 4 




Late- types 




-20.32l°j° 


-1-3018:?? 


1 7 7+5-3 
l'-'_4.4 




Starbursts 




Fit to the same faint-end slope 










All types, a = —1.37 








0.1 < z < 0.5 


91 9 9+0.10 
z±.zz_ ()6 


-, 07+0.02 

1 -°'-0.01 


28.1+?' 8 


26.73±0.04±0.06 


All 


0.5 < 2 < 0.75 


-21 7q+ - 04 




20.2l° :4 


26.81±0.03±0.07 


All 


0.75 < z < 1.0 


-21 62+ ' 04 
zr.oz_ 04 




19.5l°j 


26.73±0.05±0.07 


All 



Note. — (1) i?-band selected sample. 
(2) Luminosity densities are derived by integrating LFs characterized by the best-fitting Schechter function parameters. First 
errors represent statistical errors, while second errors represent cosmic variance (see text for details). For early- types marked 
with an asterisk (*), we have excluded the faint up-turn at Mb > —17.5. 
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Table 4. Best-fitting Schechter function parameters for rest frame J-band luminosity 

functions 



Redshift range Mj — 5 log hro 



10- 4 (Mpc//i 70 )- 3 mag" 



crg s 



- 1 Hz-^hTo/Mpc 3 ) 



Spectral types 



0.1 < z < 0.5 



0.5 < z < 0.75 
0.75 < z < 1.0 



-23.68lg;|| 
-22.97l?;£ 
-22.44l»j° 

-23.29lS:£ 

-20 71 + ' 86 
zu -' 1 -1.00 
oo oo+0.12 

-2z.m_ 016 
-23.09t°J4 



-1.48 
-0.74 
-0.45 
-1.37 
-1.65 
-1.09 
-1.31 



0.06 
0.05 
+0.29 
0.25 
+0.39 
0.33 
+ 0.08 
0.06 
+0.31 
0.24 
+0.06 
0.07 
+0.10 
0.09 



All types 
7 7+3.7 
'•'-2.8 

8.6 
12.7 

8.2 
9.4 
40.5 
19.7 



+4.6 
4.2 
+3.8 
5.1 
+4.2 
-3.0 
18.0 
7.0 
+6.5 
,.9 
+6.3 
4.9 



27.23±0.06 



27.42±0.04 
27.30±0.04 



All 

Early-types 
Early-types* 
Late- types 
Starbursts 

All 

All 



0.1 < z < 0.5 
0.5 < z < 0.75 
0.75 < z < 1.0 



-23.68tj>;£ 

oo 77+0.22 
zo. I I _ 22 

— 23 42+ - 22 



-1.48 



Fit to the same faint-end slope 

All types, a = -1.48 
+0.06 7 7+3.7 

'•'-2.8 

10-4111 

12.2±i;i 



-0.05 



27. 


23±0 


06 


27 


40±0. 


.04 


27. 


32±0. 


.05 



All 
All 
All 



Note. — (1) Ks-band selected sample. 
(2) Luminosity densities are derived by integrating LFs characterized by the best-fitting Schechter function parameters. Errors 
represent statistical errors. For early-types marked with an asterisk (*), we have excluded the faint up-turn at Mj > —19.0. 
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Table 5. Best-fitting Schechter function parameters to rest frame U-, B- and it-band 

luminosity functions 



Redshift range 


— 5 log h 7a 


a 


10 


^(Mpc/feTo)- 3 


mag 1 


erg s -1 Hz _1 (h 7 o/Mpc 3 ) 


0.10 < z < 0.82 
0.82 < z < 1.12 
1.12 < z < 1.37 
1.37 < z < 1.59 
1.59 < z < 1.80 
1.80 < z < 2.00 


-20.28+°; 4 

-2o.26t ;;° 

-20.301°;;° 
-20.48l°;i 
-20.69l°|« 
-21.16l°; 3 3 « 


-1 20+°-° 7 




54.3+1-2 
55.2l 6 5 ; i 
50.5.61^ 

9 c- o + 14.2 

18.311? 




26.55±0.03 
26.55±0.04 
26.53±0.05 
26.45±0.06 
26.38±0.13 
26.43±0.16 


Redshift range 


M* - 5 log h T0 


a 


10" 


^(Mpc/fero)- 3 


mag -1 


l°g(p^) 
erg s -1 Hz _1 (h7o/Mpc 3 ) 


0.10 < z < 0.82 
0.82 < z < 1.12 
1.12 < z < 1.37 
1.37 < z < 1.59 
1.59 < z < 1.80 
1.80 < z < 2.00 


-21 43+ u - 14 

ZL.10_ Q 13 

-21.45!°;;° 

-21.48l°;; 2 2 
-21.60l°;« 
-21.48l°i« 
-21.48l°j« 


-, 9 o+U.Ub 




37.91-: 
35.6±|;1 
28.6+lt 
21 5+ 4 ' 1 
91 1+22.6 
zl ' 1 -10.1 
9C . k+27.8 




26.89±0.04±0.07 
26.87±0.04±0.10 
26.79±0.05±0.12 
26.71±0.06±0.13 
26.66±0.09±0.15 
26.74±0.14±0.16 


Redshift range 


M* R -b log h 70 


a 


10' 


- 4 (Mpc/fe 70 )- 3 


mag -1 


l°g(p„) 
erg s- 1 Hz-^hyo/Mpc 3 ) 


0.10 < z < 0.82 
0.82 < z < 1.12 
1.12 < z < 1.37 
1.37 < z < 1.59 
1.59 < z < 1.80 
1.80 < z < 2.00 


-22.381°^ 

-22.40l° : ;g 

-22.38l° : ;* 

— 22 45+ ' 18 

-22 49+ ' 30 
zz.*y_ 34 

99 9 o+0.40 


-1 30 +U ' U4 




23.oi;;| 

18.9+^-g 
13.9+2' 4 
l0.6l 4 ;J 
13.lU° 6 3 




27.15±0.03 
27.07±0.03 
26.98±0.04 
26.87±0.05 
26.77±0.07 
26.36±0.11 



Note. — (1) ifs-band selected sample. 

(2) Faint-end slope is fixed to value derived in lowest redshift bin. 

(3) Luminosity densities are derived by integrating LFs characterized by the best-fitting Schechter function 
parameters. Errors in luminosity density represent statistical errors. For the B— band we also represent errors 
due to cosmic variance (see text for details). 



